Electric signature of magnetic domain- wall dynamics 
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Current-induced domain-wall dynamics is studied in a thin ferromagnetic nanowire. The domain-wall dynam- 
ics is described by simple equations with four parameters. We propose a procedure to unambiguously determine 
these parameters by all-electric measurements of the time-dependent voltage induced by the domain-wall mo- 
tion. We provide an analytical expression for the time variation of this voltage. Furthermore, we show that the 
measurement of the proposed effects is within reach of current experimental techniques. 
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Introduction. Recently, applications for future memory 
and logic devices, as well as important fundamental physics 
questions, have stimulated a number of experimental and 
theoretical 9-1 1 studies of the current-driven domain wall (DW) 
dynamics in ferromagnetic nanowires. It has been shown that 
DWs can be moved by a current either parallel^— or perpen- 
dicular to the wire i 7 ' 10 ' 11 In some of the experiments short 
current pulses were employed to depin a DW from pinning 
sitesj 2 ^ Furthermore, the topological electromotive force in- 
duced by DW dynamics in a vortex DW has been studied both 
experimentally and theoreticallyi&i 2 . 

A conventional experimental method to study the DW dy- 
namics in nanowires is to measure the average DW veloc- 
ity using Kerr polarimetry, 13 x-ray microscopy,— or electron 
microscopy. 514 These types of experiments require a compli- 
cated setup which is separate from the one needed for the DW 
manipulation. This situation is neither ideal for studies of DW 
dynamics nor for further technological advances. 

In this Letter we propose a way to use the same experimen- 
tal setup for both current DW manipulation and simultaneous 
measurements of DW dynamics. Our main results are that the 




time-dependent voltage induced by the DW motion J£ can 
be used to fully and comprehensively determine the effective 
parameters of the DW dynamics. This proposal follows from 
the fundamental properties of the current-induced DW mo- 
tion, namely: (i) Applied DC current (above critical value) 
produces voltage with AC components, (ii) Applied AC cur- 
rent induces phase shifted AC voltage. The magnitude of the 
proposed effects is calculated to be within current experimen- 
tal resolution. 

Similar techniques have already shown promise in magnetic 
field driven DW systems. 17 This method should make it more 
feasible to utilize DW dynamics for device applications. Fur- 
thermore, the proposed systematic approach can be used to 
compare the extracted phenomenological parameters of the 
DW dynamics for a system described by arbitrary underlying 
Hamiltonian to those of microscopic theories. 

Model. The dynamics of the magnetization S in a quasi- 
one-dimensional wire is described by Landau-Lifshitz-Gilbert 
(LLG) equation with current jr^- 

S = -S x H e -jd z S + (3jS x d z S + aS x S, (1) 

where H e = —6H/SS is the effective magnetic field given 
by the Hamiltonian H of the system, S = M/|M| is a unit 



FIG. 1 . (Color online) A moving head-to-head domain wall of width 
A. The DW is centered at zq and is tilted by an angle <j>. 



magnetization vector, a is the Gilbert damping constant, j3 is 
the non-adiabatic spin torque constant, d z = d/dz where z 
is along the wire, and the time is measured in units of the gy- 
romagnetic ratio 70 = g\e\/ (2mc). DWs in a ferromagnetic 
wire can be modeled by a spin Hamiltonian H which contains 
exchange, spin-orbit, 20 and dipolar interactions. In a thin wire, 
the latter can be approximated by two anisotropies: a strong 
anisotropy along the wire (A) and a weak anisotropy trans- 
verse to it (K). In realistic systems a, j3 -C 1 and K <C A. 

In a thin wire, a lowest-energy magnetization configuration 
(at j — 0) is uniformly ordered along the z or — z direc- 
tion. A static DW is the next low-energy configuration with 
the boundary conditions 5 z (±oo) = ±1 or S^ioo) = =pl, 
DWs can be injected in the wire using different techniques. A 
sketch of a wire with a DW of width A, determined by the 
Hamiltonian parameters, is depicted in Fig.Q] 

For small enough applied currents, it can be shown that the 
DW in a thin wire is a rigid spin texture 14 and its dynamics can 
be described in terms of only two collective coordinates.^^ 2 . 
These coordinates correspond to the two softest modes of the 
DW motion: the DW position along the wire, Zq, and the rota- 
tion angle <j> of the magnetization in the DW around the wire 
axis, see Fig. Q] It has been shown 22 - 23 that the equations of 
motion for the DW in a thin ferromagnetic wire are model 
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independent and can very generally be written in the form 

Zo=Aj + B[j-j c sm(2<f>)} ) (2) 
j> = C\j-j c sin(2<f>)]. (3) 

Here all current nonlinearities are neglected, since the large 
currents leading to observable nonlinear effects would burn 
the nanowire. For a dc current below the critical value j c , i.e., 
j < j c , Eq. (01 implies that the DW tilts from the transverse 
anisotropy plane by the angle that satisfies sin(2(/>) = j/j c 
around the wire axis and then moves along the wire with a 
constant velocity Aj. For j > j c , the DW constantly rotates 
while moving. 

The coefficients A, B, C and the critical current j c are 
the parameters that fully describe the DW dynamics. They 
can be calculated microscopically for certain toy models, 22 
but in general they vary for different wires and depend on 
the temperature and nanofabrication details. Therefore, in 
this Letter we propose a way to determine these coefficients 
by model-independent measurements of an induced ac volt- 
age directly from an experiment suitable for all-electric DW 
manipulation. As we show below, this ac voltage can be in- 
duced by applied dc currents and by certain time-dependent 
current pulses with parameters similar to those achieved in re- 
cent experimentsi^i 2 ^ 

Microscopically the dynamics parameters can be obtained 
in the following way. The energy of a static DW, Eq(zq, qb) — 
J H[Sq(z, zo, <p))dz, where So is a solution of a static LLG 
with K = 0, in general depends on both zq and <fi. However, 
assuming that the wire is translationally invariant (pinning can 
be neglected), Eq would not depend on the DW position zq 
and therefore d Z(l Eo — 0. The only contribution to Eq that 
depends on the angle <j> comes from the small anisotropy in 
the transverse plane, Er,(<f)) = — k cos(20) ( 2 2 ^ This allows 
us to find the coefficients in Eqs. (O and (f3]l in terms of the 
parameters of the LLG (Q3'^ i ^ Up to first order in a and B 
they are 

B „ a — 6 , 
A=^, B = —r^{l + aa z4> ) 7 (4) 
a a 
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where a = aV, ft = 0T>, V = 

a zz = I / dz(d z S Q ) 2 , a u = \ J dz^So) 2 , and a z<t> = 
\ J dzd z So-d<f,So. Equations (HJi and (O are consistent 2 ^ with 
the expressions for A, B, C, and j c found in Ref.l22i 

We now outline the method to find A, B, C, and j c directly 
from all-electric measurements. It is based on measuring the 
ac voltage V induced by a moving DW. To find V one has to 
know the time evolution of the total energy (per unit area of 
the wire's cross-section) in the system, 
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In general, DW energy has two contributions: the power sup- 
plied by an electric current and a negative contribution due to 
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FIG. 2. (Color online) Dependence of average voltage (V) on dc 
current j for C > and C < 0, respectively, see Eq. ([§}. The slope 
at j < j c gives ^g^, whereas at j S> j c it gives + (1 + A)C. 



dissipation in the wire. Using the general solution of the LLG, 
Eq. (HJ, one can obtain the derivative of the energy a a 22 ' 27 



E = 2[0a zz z a + (1 - j3a z(t> )(j)\j - a 



dzSg. 



(7) 



The last term on the right-hand side of Eq. ((7) describes the 
dissipation and is therefore always nonpositive. Meanwhile, 
the first term is proportional to the current density j and gives 
the power Vj supplied by the current. With the help of 
Eqs. ©-(O and adopting the approximation V ~ 1 of Ref.l22l 
we obtain the expression for the induced DW voltage 29 , 
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j + C(l+A)\j-j c sm(2<j>)]. 



(8) 



Note that Eq. ([8]) gives the contribution to the voltage due 
to DW motion. This contribution is in addition to the usual 
Ohmic one. The voltage V in Eq. © is measured in units 
of PgnB /(e7o) and the current density is measured in units 
2eMj (Pg[iB), where P is the current polarization. We em- 
phasize that unlike in the previously studied cases,^^ this 
voltage is not caused by the motion of topological defects 
(vortices) transverse to the wire. 

Measurement of coefficients A, B, C, and j c . In order to 
find coefficients A, B, and C, we propose three independent 
measurements of the voltage induced by a moving DW. Al- 
though there are various factors affecting the nanowire resis- 
tance, the contributions from most of them are independent 
of DW motion and therefore give only a constant component 
of the resistance. To characterize the DW dynamics, one has 
to concentrate only on the resistance variations in time. Our 
estimates show that the amplitude of voltage oscillations due 
to DW motion is of the order of 10~ 7 V and therefore experi- 
mentally measurable. 

Equation ([8]) implies that the voltage of the DW can give all 
the necessary information about DW dynamics. Namely, one 
can obtain C by measuring the voltage changing with time 
and parameters A and B by measuring the amplitude of the 
voltage oscillations. 

Slopes measurement. In Refs.l23landl30 it was proposed to 
obtain A, B, and j c by measuring the drift velocity of the DW, 
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(io). It is important to note that Eq. (0 has the same form as 
Eq. (f2j). Thus, instead of measuring the drift velocity, which 
requires a more complicated experimental setup, we propose 
to perform all-electric measurements. Namely, to measure the 
average voltage of DW, (V), as a function of dc current. From 
Eq. © one can see that (V) = ^jp-j for j < j c , whereas 

(V) = ^rJ + il + A^^f - ]? for j > j c , seeFig.E] The 
critical current is determined by the end of the region linear in 
j for small currents. The measurement of slope k\ at j < j c , 
and slope k 2 at j ^> j c gives the two independent quantities: 



V(f)[ 



C* > 



C < 
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k 2 -*! = (! + A)C. 



(9) 



Instead of measuring voltage average for dc current, one can 
apply a linearly increasing time-dependent current j(t) = qt 
below the critical value j c . At sufficiently small q the voltage 
will also be linear in time, V(t) w ^jf-qb. By measuring this 
voltage one can find 



V(t) _ tfC_ 



B 



(10) 



Once C is determined, Eqs. © give A and B. The drawback 
of this measurement is that it might be hard to disentangle ki 
and k 2 from the Ohmic contribution. However k 2 — fci is free 
from the Ohmic resistance of the wire. 

In order to find C, the most intuitive approach is to input a 
dc current slightly above j c . Then the voltage induced by the 
moving DW will oscillate with the period of the double angle 
</>, see the insets of Fig. [3] The half- width of the peak (dip ) 
for C > (C < 0) is given by arccos(j c /j)/(|CVj 2 - ft). 
The measurement of the voltage oscillations period Tq (which 
we estimate to be ~ 10~ 7 - 10 -6 s) determines C at a given 
j- 



I JQ 



j-j c sin(20) T ^f^l 



(11) 



F° r j ~ jc "C j c > the period diverges but the half-width 
~ l/(Cj c ) stays finite. To obtain the period To, one can 
perform the Fourier transform of V(t) to find the frequency 
/o = l/T Q , seeFig.E] 

To determine coefficient A in the same experiment, one can 
measure AV — V max — V m i n = 2(1 + A)|C|j c , see insets of 
Fig. [3] Then 



A = 



AV 

2\C\Ic 



- 1. 



(12) 



Note that AV = 2(k 2 — k\)j c and therefore this experiment 
can also provide a crosscheck with the aforementioned mea- 
surement of the slopes. 

Phase shift experiment. Another method to measure the 
coefficient C is by applying an ac current j = jo sin uit with 
jo > jc, which has only a short time interval where j > j c , 
so that there is only one period of voltage within the period 
of j(t). One can measure the phase delay, AO, between the 
current maximum and voltage extremum 31 (see Fig. [4]). Next, 
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FIG. 3. (Color online) Fourier transform of the voltage V as a func- 
tion of frequency / at the dc current 1.1 j c . The insets show V as a 
function of time t for C > given by a = 0.02 and /3 — 0.01; and 
for C < given by a = 0.01 and j3 — 0.02. The voltage period 
is To = 1/fo- In the inset for C < 0, the voltage varies between 
V max = 0.041j c /A and V mm = -0.019j c /A. 





(a) 




v max 








/ 

/ 




\ \ 
\ \ 




/ 

/ 




V \ 





-» A0 






( c ) 


- - - j 




y 


x 

V 




/ ■ ■ 

S 

/ 

/ 

/ 




\ 

\ 

\ 










v;„„ 



(b) 






/ 

/ 

/ ^. 




\ 

\ 




A0 


(d) 








/ 

/ 






\ 







U)t 



FIG. 4. (Color online) Input current j (dashed line) and measured 
voltage V (solid line) as functions of time t. (a) and (b) show 
the phase delay A# between the current maximum and voltage ex- 
tremum for C > and C < 0, respectively, (c) and (d) depict V(t) 
at A9 — for the same C > and C < 0, respectively. 



one fixes the amplitude jo and tunes the frequency ui until 
A8 = 0. In this case, for jo — j c <C j c , we can use half 
of the time interval for which the current pulse is above j c to 
approximate the period of <p by dc current j as 



1 I 9 ' 3° 

— 7r — 2 arcsm — 
2w V jo 



IcIx/jo^jF 



(13) 



For jo — j c <C j c , Eq. ( fT3l can be further simplified to give 

\C\ « (14) 
2(jo - Jc) 

In other words, when w w C(jo — j c ) which corresponds 
roughly tow~ 10 7 Hz, the current pulse covers only one pe- 
riod of voltage. Our simulations show that the expression (1141 ) 
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see Fig. [5] To measure the decay of V (i) with time, one inputs 
a dc current below j c , then measures voltage V] immediately 
after turning off the current at ti, and then later measures volt- 
age Vf at time i; + A t . We note that right after turning off 
the current, there is a short time period when the DW dynam- 
ics cannot be described by Eqs. <J2J and (0. It corresponds to 
the dynamics of fast degrees of freedom. This process has a 
characteristic time ~ 10~ n s which is typically much smaller 
than the voltage decay time ~ 10~ 8 s. Thus we can safely 
assume that the rotation angle cf) does not change much during 
this time interval, and we find 



|C| 



1 



■In 



2Vi/V f 



2A tJc 1 + y/l-p/jl 



(15) 



FIG. 5. (Color online) Voltage (solid line) evolution after the current 
(dashed line) is turned off at time ti for C > given by a — 0.02 
and ft — 0.01. Inset: the same dependencies for C < given by 
a = 0.01 and /3 = 0.02. The measurement of Vf is performed at 
ti + At. The region encircled by the dotted line cannot be described 
within our approach but it is too small to effect our results. 

works sufficiently well for jo < 1.3j c . The sign of C is deter- 
mined by the extremum of the measured voltage: C > if V 
has the minimum and C < if V has the maximum. 

Our simulations show (Fig. @|l that in addition to the large 
peak (dip) of voltage there is a smaller one with the opposite 
curvature. This is because when j(t) reaches j c , the angle cf> 
has not yet rotated to the angle corresponding to sin(2<^rj) = 1 
due to the cumulative phase delay between current and volt- 
age. 

Abrupt current pulse experiment. It is also possible to mea- 
sure the coefficient C for currents below the critical value j c . 
The constant \C\j c determines the internal time scale of the 
DW motion. After one switches the subcritical current off at 
time ti, the voltage asymptotically decays as exp(— 2\C\j c t), 



which is valid for Vi/Vf 1. For example, estimating 
Vi/Vf = 10 we find |C| « 1.17/(A t j c ). The sign of C can 
be easily determined by the form of voltage decay (see Fig. |5j. 

To summarize, we propose several all-electric measure- 
ments of the parameters fully describing domain-wall dynam- 
ics in thin ferromagnetic nanowires. These measurements are 
based on the voltage induced by a moving DW in response to 
certain current pulses. Our proposal opens doors for experi- 
ments which are suitable not only for all-electric DW manip- 
ulation but also for the simultaneous measurement of the DW 
dynamics. These findings give a more reliable and straight- 
forward experimental method to determine the DW dynamics 
parameters, which can then be compared to microscopic the- 
ories. The procedure we described works for a given temper- 
ature regime. It may also be used to investigate the tempera- 
ture dependence of the effective parameters. Future work will 
include accounting for pinning effects, which brake transla- 
tional invariance in the wires. 32 
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